Isospin dependence of pseudospin symmetry in nuclear resonant states 
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The relativistic mean field theory in combination with the analytic continuation in the coupling 
constant method is used to determine the energies and widths of single-particle resonant states in 
Sn isotopes. It is shown that there exists clear shell structure in the resonant levels as appearing 
in the bound levels. In particular, the isospin dependence of pseudospin symmetry is clearly shown 
in the resonant states, is consistent with that in the bound states, where the splittings of energies 
and widths between pseudospin doublets are found in correlation with the quantum numbers of 
single-particle states, as well as the nuclear mass number. The similar phenomenon also emerges in 
the spin partners. 

PACS numbers: 21.10.Hw, 21.60.Jz, 25.70.Ef, 27.80.-l-w 



INTRODUCTION 

About 30 years ago a quasi degeneracy was observed in 
heavy nuclei between single-nucleon doublets with quan- 
tum numbers (n, l,j = 1+1/2) and {n—1, 1+2, j = 1+3/2) 
where n, I, and j are the radial, the orbital, and the total 
angular momentum quantum numbers, respectively 
The quasi-degenerative states are suggested as the pseu- 
dospin doublets j = I ± s with the pseudo orbital an- 
gular momentum I, and the pseudospin angular mo- 
mentum s, and have explained a number of phenom- 
ena in nuclear structure including deformation 3] and 
superdeformation0 , magnetic moment 0, Q and iden- 
tical bands0, 0, IJ. Because of these successes, there 
have been comprehensive efforts to understand its ori- 
gin since the discovery of this symmetry. Blokhin et 
al.T^ performed a helicity unitary transformation in a 
nonrclativistic single-particle Hamiltonian, showed that 
the transformed radial wave functions have a different 
asymptotic behavior, implying that the helicity trans- 
formed mean field acquires a more diffuse surface. Ap- 
plication of the helicity operator to the nonrclativistic 
single-particle wave function maps the normal state {I, 
s) onto the pseudo state {I, s), while keeping all other 
global symmetries. The same kind of unitary transfor- 
mation was also considered earlier by Bahri et al. to 
discuss the pseudospin symmetry in the nonrclativistic 
harmonic oscillator They showed that a particular 
condition between the coefficients of spin-orbit and orbit- 
orbit terms, required to have a pseudospin symmetry 
in that nonrclativistic single-particle Hamiltonian, was 
consistent with relativistic mean-field (RMF) estimates. 
Based on the RMF theories, Ginocchio identifies a possi- 
ble reason for this; namely that the symmetry arises from 
the near equality in magnitude of an attractive scalar, S, 
and repulsive vector, V, relativistic mean field, S ~ —V, 
in which the nucleons move[l^. They reveal the pseu- 
dospin symmetry is exact when doublets are degener- 
ate. Meng et al.yjjllJI show that pseudospin symmetry 
is exact when c?S/dr — where 'S = V + S and the 
quality of the pseudospin approximation in real nuclei 



is connected with the competition between the pseudo- 
centrifugal barrier and the pseudospin-orbital potential. 
In other works 0, 0, the pseudospin symme- 
try has also been discussed in relation with the argu- 
ments and conjectures of Refs.fl^ H^ . Pseudospin 
symmetry has also been shown to be approximately con- 
served in medium energy nucleon scattering from even- 
even nuclei 0, 0|. A test of nuclear wave functions 
for pseudospin symmetry was done in Ref . |20l [2l| . The 
structure of radial nodes occurring in pseudospin levels 
and a classification for the intruder levels states which do 
not have a pseudospin partner in the limit of pseudospin 
symmetry was explained as a direct effect of the behav- 
ior of nodes of Dirac bound states |23|, giving support for 
the relativistic interpretation of nuclear pseudospin sym- 
metry. Pseudospin symmetry is also investigated for the 
relativistic harmonic oscillator, and relativistic Woods- 
Saxon and Hulthen problem, in which the analytical solu- 
tions are obtained for the bound states of the correspond- 
ing Dirac equations by setting either T, or A {= V — S) 
to zerolH lil m m m ill. The pseudospin breaking 
is shown to be connected with the mean-field shaped by 
the parameters of the harmonic oscillator potential[22|, 
or the Woods-Saxon potentialjlSl- All these indicate 
the pseudospin symmetry is a good approximate for the 
single-particle bound states of atomic nuclei. However 
for the resonant states, which are thought to play a im- 
portant role in formation of the exotic phenomena such 
as halo and skin[3lL 13^. l33j|. the symmetry has not been 
examined in detail for real nuclei. In recent years, many 
techniques have been developed to study resonant states, 
such as the R-matrix theory jSJLthe extended R-matrix 
theory [23, the K-matrix theor y l3fil |. S-matrix theory j^^j 
the real stabilization method j38i|. the complex scaling 
methodli^j and so on. Compared with these meth- 
ods, the analytic continuation in the coupling constant 
(ACCC) method proposed by Kukulin[4o| is more effec- 
tive and numerically quite simple, has been applied to 
investigate the energies and widths for resonant states 
in the framework of the non- relativistic 113, E3| and 
relativistic theory[3E3- Based on the Schrodinger and 
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Dirac equations, we have investigated the stabihty and 
convergence of the energies and widths for single-particle 
resonant states with square-well, harmonic oscillator and 
Woods-Saxon potentials, and discussed their dependence 
on the coupling constant interval and the Pade approx- 
imant (PA) order 0, ji^l- Particularly, we have stud- 
ied the pseudospin symmetry in the resonant states for 
Dirac- Woods-Saxon problem, and shown the correlations 
of pseudos pin splittings with the shapes of Woods-Saxon 
mean field '4q|. Considering that the RMF theory is a 
more microscopic one, and has gained considerable suc- 
cess in describing many nuclear phenomena for the stable 
nuclei as well as nuclei even far from stability, here we 
intend to study the resonant states including their pseu- 
dospin and spin symmetries by the ACCC method in the 
framework of RMF theory, as done for the bound states 
in the RMF theory 14]. 



THEORETICAL FRAMEWORK 

The basic ansatz of the RMF theory is a Lagrangian 
density whereby nucleons are described as Dirac parti- 
cles which interact via the exchange of various mesons 
(the scalar cr, vector uj and iso- vector vector p) and the 
photon ^3 

C = - M)^P + Id^ad^'a - U{(j) - \VL^,^VL^''' 
+ \mlLO^LO^^ - iR^.R'^'' + im^p^pM _ iF^^Ff^- (1) 
-ga'ipcril; - guo^^tpil) - gpi'iJTip - e-0^V', 

where M is the nucleon mass and ma {go), rriui (ffw), 
and nip (gp) are the masses (coupling constants) of the 
respective mesons. A nonlinear scalar self-interaction 

has been included [5C 
mesons are given as 



-cr^ of the a meson 



The field tensors for the vector 
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df^p" - d^pf" ~ gPipf" X p") 



(2) 



The classical variation principle gives the following 
equations of motion 



[a-p + W(r)+/3(A/ + l^s(r))]^ = eV-, 



(3) 



for the nucleon spinors, where e and -0 are the single- 
particle energy and spinor wave function respectively, 
and the Klein-Gordon equations 



(-Ao- + U'{a)) = g^ps, 

(-A + m2)co^ = .g^j^(r), 

(-A + m^jpA' ^ .gpj^(r), 

-A <(r) = ej;;(r). 



(4) 



for the mesons, where 

Vv(r) = g^ip + gp^T + ]-e{l~T3,)4^, 



Vs{v) 



(5) 



are the vector and scalar potentials respectively and the 
source terms for the mesons are 
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(6) 



It should be noted that the contribution of negative en- 
ergy states are neglected, i.e., the vacuum is not polar- 
ized. Moreover, the mean field approximation is carried 
out via replacing meson field operators in Eq. (PJ by 
their expectation values, since the coupled equations Eq. 
Q and Eq. Q are nonlinear quantum field equations 
and their exact solutions are very complicated. In this 
way, the nucleons are assumed to move independently in 
the classical meson fields. The coupled equations can be 
solved self-consistently by iteration. 

For spherical nuclei, the potential of the nucleon and 
the sources of meson fields depend only on the radial 
coordinate r. The spinor is characterized by the angular 

momentum quantum numbers I, j, to, the isospin t — ± — 

for neutron and proton respectively. The Dirac spinor 
has the form 



V'(r) 



Xtit), 



(7) 



where Yj^{6, (/>) are the spinor spherical harmonics. The 
radial equation of the spinor, i.e. Eq. |(2Jl, can be reduced 

to 113 

(-|: + ^)Fi, (r) + [M + Vpir)] (r) = eGy (r), 

(|: + (0 - [M ~ V^{r)] Fi, (r) = eFi, (r), (8) 

in which Vp{r) = Vv{r) + Vs{r), Kn(r) = W(r) - Vs{r), 
and K = (— l)-'+'+^/^(j' + 1/2). The meson field equations 
can be reduced to 



Qj.2 



2d_ 

r dr 



+ m! U) = s^(r). 



(9) 



where = cr, w, p, and photon ( = for photon) 
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The source terms read 

520'^ (r) — .93(7^ (r) for the a field, 



9ujPv 



for the Lu field, 
for the p field, 
for the Coulomb field, 



In the non-relativistic limit (fc ^ M), Eq. H14|l reduces 
to 

4-2 _ 1,2 Oh h 

(15) 



E = 



(10) 



2M 



r = 



M 



and 
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n=l 

i:(|G,(r)|2 + |F,(r)n. 



(11) 



By solving Eqs. (|HJ) and (|5Jl in a meshed box of size 
Ro self-consistently, one can calculate the ground state 
properties of a nucleus. The vector potential Vv{r) and 
the scalar potential Vs{r), energies and wave functions 
for bound states are also obtained. By increasing the 
attractive potential as Vp{r) — > AV^(r), a resonant state 
will be lowered and becomes a bound state if the coupling 
constant A is large enough. Near the branch point Aq, 
defined by the scattering threshold fc(Ao) = the 
wave number A: (A) behaves as 



fc(A) 



iVA - Ao, I > 0, 
i(A-Ao), ^ = 0. 



(12) 



These properties suggest an analytic continuation of the 
wave number k in the complex A plane from the bound- 
state region into the resonance region by Fade approxi- 
mant of the second kind (PAII) 40] 



k{x) w k^^'^\x) 



Co + CiX + C2X + . . . + ClX 



1 + dix - 



(13) 

where x = \/\— Aq, and cq, ci, . . . , cl, di, (i2, . . ■ , c^Af are 
the coefficients of PA. These coefficients can be deter- 
mined by a set of reference points Xi and k{xi) obtained 
from the Dirac equation with Ai > Aq, i ~ 1, 2, L-f A'^-f 
1. With the complex wave number fc(A — \) = k^ + ik^, 
the resonance energy E and the width F can be ex- 

F 

tracted from the relation e = E — i— (E,T e 



and 



I.e. 



E 



^{KP + fc2 - fc2)2 + 4fc2fc2 ^ (^2 ^ ^2 _ ^2) 
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Fig. 1: The single-particle energies of neutron as the function 
of mass number for Sn isotopes, where the energies of resonant 
states obtained by the ACCC method in the RMF theory with 
the interactions NL3. 
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Fig. 2: The same as FigQ but for the single-particle widths 
of neutron. 
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NUMERICAL DETAILS AND RESULTS 



With the technique represented above, we first single 
out the resonant states in the continuum for Sn isotopes, 
and then analyze their properties including the pseu- 
dospin and spin symmetries. Following Refs. [isl lia. lirl 
li^ . the radial Dirac equations © are solved in a stan- 
dard way by the shooting method with a step size of 0.05 
fm using proper boundary conditions in a spherical box of 
radius R — 20 fmji^. The branch point Xq is determined 
in a self-consistent way as in Ref. |40j] . In correspondence 
to L + N + 1 different coupling constant values Ai > Ao, 
we solve the bound state problem for Ail](r), and ob- 
tain the eigenvalues Si from the Dirac equations. The 
solutions of the linear equations are the coefficients in 
the PA expression With the coefficients in Eq.lUSJ, 
the wave number as a function of the coupling constant 
fc(A) can be obtained. Then the energy and width of a 
resonant state can be extracted from the complex wave 
number k{X — 1). Considering that there are a little in- 
fluences of the coupling constant interval and the Fade 
approximant order on the calculated results in the ACCC 
method as shown in Ref.ji^, we first check the stability 
and convergence of this method in the present framework 
to choose reasonably the coupling constant interval and 
the Fade approximant order for Sn isotopes. Similar to 
the Figs. 2 and 3 in Ref. [isj . the energies and widths of 
resonant states are found to be considerable stable in a 
fairly large coupling constant interval for fixed order, es- 
pecially for the calculations with (5,5) FAs. Based on the 
fact, all the energies and widths for resonant states in the 
following are taken from the calculations of (5,5) FA with 
Af, — Ao = 1 as done in Ref. [48| . in which the branch point 
Ao and Ab are respectively the left and right edge of the 
coupling constant interval in the Fade approximant. 

To see the behavior of the single-particle resonant 
states and their isospin dependence as well as compar- 
ison with the case of bound states, the energy as a func- 
tion of the mass number are displayed in Fig.l for the 
Sn isotopes. Except for several small kinks, the energies 
are shown to decrease monotonously with the increasing 
isospin for all the resonant states, which is in agreement 
with that for the bound states near the threshold, differ- 
ent from that for the bound states far from the threshold 
where the change is altered with a contrary trend as the 
level lsi/2- For the resonant states 8^3/2,3^5/2, 8/5/2, 
and 3/7/2, the energies decrease more quickly, and hap- 
pen to cross with their neighbor at several Sn nuclei, 
which does not appear in the bound states. Particularly, 
there are several kinks in these levels as appearing in the 
bound levels at a less extent, e.g.: the lsi/2,2si/2, 2pi/2, 
and 2p3/2. These kinks always appear at several spe- 
cial positions for whether the bound or resonant states, 
may reflect some special structure existing in these nu- 
clei. Besides of the kinks, several large gaps are found in 
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Fig. 3: The energy splittings of pseudospin doublets as the 
function of mass number for Sn isotopes, where the filled and 
opened marks represent respectively for the cases of bound 
and resonant states. 



25 H-' ' I 



I ' ' ' ' I ' ' ' ' I ' ' ' ' I- 

<^-E(3d5/2)-E(2g7/2) 
■i— E(3f7/2)-E(2h9/2) 
-<^E(2g9/2)-E(l ill/2) ■ 




160 



170 



110 120 130 140 150 
Mass number 

Fig. 4: The same as Fig. 3, but for the width splittings of 
pseudospin doublets for Sn isotopes. 



the resonant levels as appearing in the bound levels. The 
large energy gap between the levels 2/7/2 and I/111/2 re- 
main hardly change whether these levels are unbound or 
lowered to bound. The same phenomenon also appear in 
the resonant levels between lji3/2 and 2/ig/2. From these 
phenomena we see that the resonant levels also provide 
the detailed information on nuclear shell structure. 

Compared with the energies, the widths of the reso- 
nant states also decrease monotonously with the increas- 
ing isospin, which can be seen in Fig. 2. Especially for 
the resonant states 8^3/2, 8^5/2, 8/5/2, and 8/7/2, the 
behavior of decreasing more quickly, even crossing with 
their neighbor in the energies also appears in the widths. 
Moreover, the shell structure is clearly shown with sev- 
eral large gaps in the widths for several different reso- 



5 



nant states. The large energy gap between 2/ig/2 and 
3/7/2 decreases with the increasing isospin. The same 
phenomenon also appears between the levels liii/2 and 
259/2- In comparison with the energies, the widths de- 
crease in general more quickly, especially for the 8/5/2 
and 3/7/2. 

With these knowledge on the energies and widths, we 
analyze the symmetry existing in these resonant states. 
As identified to be a good approximation in the bound 
states, the pseudospin symmetry worths being examined 
for the single-particle resonant states. In R,ef . [i^ . the 
symmetry in the resonant states has been investigated 
by solving the Dirac equation with Woods-Saxon po- 
tentials in combination with the ACCC method, where 
the pseudospin breaking is shown in correlation with the 
shapes of nuclear mean-field. In order to examine the 
symmetry in the framework of RMF theory, and study 
their isospin dependence, which has been discussed in 
the bound statesjQ, we analyze the relation between 
the pseudospin splitting and the nuclear mass number for 
the resonant states. The dependence is shown in Fig. 3 
for the Sn isotopes, where the energy splittings of bound 
states are also exhibited for comparison. From Fig. 3, it 
can be seen that the pseudospin symmetry in the reso- 
nant states is correlated with the quantum numbers of 
single-particle states, and the pseudospin splittings are 
different from the different pseudospin partners, which is 
in agreement with the case of bound states [lOI as well as 
the case of resonant states in Dirac Woods-Saxon mean 
field approximation^^. As far as the isospin dependence 
of the pseudospin orbital splitting is concerned, the split- 
tings between the resonant doublets give a monotonous 
decreasing behavior with the increasing isospin with only 
the 3/7/3 and 2/ig/2 partner exception. Particularly for 
the 2(7g/2 and liii/2 partner, the pseudospin splitting in 
-^^^Sn is only half of that in -^°^Sn. Just as we expected, 
the pseudospin symmetry in neutron-rich nuclei is bet- 
ter, which agrees the case of the bound levels. Besides of 
the isospin dependence, the pseudospin splitting is also 
correlated with the radial, orbital, and total angler men- 
turn quantum numbers. Such as the splitting between 
the 2/111/2 and IJ13/2 partner is 2 MeV larger than that 
between the 2179/2 and liii/2 partner. In comparison 
with the resonant states, the situation in the bound pseu- 
dospin partners is more complicated, but the pattern is 
more or less the same, i.e., a decreasing behavior with 
the increasing isospin with some exceptions. For exam- 
ple, the splittings between the levels 2sx/2 and 1^3/2, 
and between the levels 2^3/2 and I/5/2 show a decreas- 
ing trend in a complicated way, while that 2^5/2 and 
1(77/2 presents a minor variation in the whole Sn isotope 
chain. It should be mentioned for the 2/7/2 and lft.9/2 
partner, the same trend for the change of pseudospin 
splitting is seen whether they are bound or evolved to 
unbound states. For all the pseudospin partners with the 
same radial quantum number, the energy splitting is dis- 



covered to evolve from £^„,/j=/+i/2 > £'„_i,/+2,:,=i+3/2 
to -E„_/j=/+i/2 < £^„-i,;+2j=;+3/2 with the increasing 
orbital angular momentum. This inversion of pseu- 
dospin splittin gs i s observed experimentally, was found 
in Refs.[lJ,|l^|5lllS3l i'^ bound states, also appears 
in the resonant states here. 

For width (to see Fig. 4), similar to the energy, the 
pseudospin splittings decrease monotonously with the in- 
creasing mass number with no exception. The widths 
for all the single particle states with quantum numbers 
{n,l,j = 1 + 1/2) are systemically larger than that of 
their pseudospin partners ( n — 1, ^ -I- 2, j = / + 3/2), 
which implies that the resonant states with higher or- 
bital angular momentum hold longer decay time than 
their pseudospin partners. It is because that the orbital 
angular momentum is smaller, the centrifugal barrier is 
lower, the resonant width is larger. Hence, even the en- 
ergy is fully degenerate for the pseudospin doublet, their 
decay life maybe still different, which is noticeable to ex- 
plore the resonant state from experiment. From these 
studies, we see that the pseudospin symmetry remains a 
good approximation for both stable and exotic nuclei. A 
better pseudospin symmetry can be expected for the or- 
bital near the threshold, particularly for nuclei near the 
particle drip line. 
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Fig. 5: The energy splittings of spin doublets as the function 
of mass number for Sn isotopes, where the filled and opened 
marks represent respectively to the cases of bound and reso- 
nant states. 

In order to compare with the pseudospin symmetry, 
the spin orbital splittings as a function of the mass num- 
ber are plotted in Fig. 5 for the Sn isotopes, where a 
more distinct decreasing behavior with the increasing 
mass number is clearly seen for both the bound and un- 
bound spin partners with few exceptions. It shows that 
the spin symmetry in neutron-rich nuclei is also bet- 
ter such as the pseudospin symmetry. In addition, the 
spin orbital splitting is also correlated with the quan- 
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Fig. 6: The same as Fig. 5, but for the width splittings of spin 
doublets for Sn isotopes. 



turn numbers of single-particle states. For the spin dou- 
blets with the same radial quantum number, the energy 
splitting increases with the increasing orbital angular mo- 
mentum for the bound states, which is contrary to the 
resonant states. Compared with pseudospin symmetry, 
the inversion of spin splitting has not been found for 
whether the bound or unbound spin partners. The ener- 
gies for the single particle states with quantum numbers 
{n,l,j — I — 1/2) are systemically larger than that of 
their spin partners {n,l,j — I + 1/2), which indicates 
that the single-particle states with higher total angular 
momentum are more stable than their spin partners. It 
agrees the experimental facts as observed in the nuclear 
bound states. For the width, similar to the pseudospin 
symmetry, a monotonous decreasing behavior with the 
increasing isospin is seen in Fig. 6 with a bit special in 
several nuclei. The widths for the single particle states 
with quantum numbers (n, l,j = 1—1/2) are systemically 
larger than that of their spin partners ( n,l,j = 1 + 1/2), 
which implies that the resonant states with higher to- 
tal angular momentum hold longer decay time than their 
spin partners, is similar to the case of pseudospin sym- 
metry. From these results we see that the spin orbital 
splitting is also very sensitive to the quantum numbers 
of single-particle states. The spin symmetry remains a 
better approximation for the bound states with lower or- 
bital angular momentum, while for resonant states with 
higher orbital angular momentum, particularly for nuclei 
near the particle drip line. 



CONCLUSION 

The relativistic mean field theory in combination 
with the analytic continuation in the coupling constant 
method is adopted to determine the energies and widths 



of single-particle resonant states for the Sn isotopes, 
where the energies are found to decrease monotonously 
with the increasing isospin with several small kinks, is in 
agreement with the case of bound states near the thresh- 
old. In addition, the shell structure is seen clearly in 
the resonant levels with several large gaps as appearing 
in the bound levels. In particular, the pseudospin and 
spin symmetries in the resonant states are investigated 
and compared with that in the bound states, where the 
symmetries are found to be correlated with the quantum 
numbers of single-particle states, which is similar to the 
case of bound states. As far as the isospin dependence 
is concerned, the pseudospin and spin splittings between 
the resonant doublets give a monotonous decreasing be- 
havior with the increasing isospin with few exceptions. 
The same phenomenon also appear between the bound 
doublets. The trend of energy splittings with the orbital 
angular momentum is found to be consistent with the 
case of bound for all the pseudospin partners, which is 
not consistent with the spin partners. The inversion of 
energy-splitting appears in several pseudospin partners 
of resonant states, which agrees with the experimental 
observation, as well as the case of bound states and the 
nonrelativistic prediction, while the inversion does not 
appear in the spin doublets. The tendency of width split- 
tings with the parameters is displayed to be agreeable for 
all the pseudospin partners, as well as the spin partners. 
The widths for the single particle states with higher or- 
bital (total) angular momentum are systemically larger 
than that of their pseudospin (spin) partners, which show 
that even the fully degenerate pseudospin (spin) doublets 
may exist different decay time, making them worthy of 
experimental attention. 
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